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The bound states of and ^He have been calculated using the Argonne uis plus the Urbana three- 
nucleon potential. The isospin T = 3/2 state have been included in the calculations as well as the 
n-p mass difference. The ^H-®He mass difference has been evaluated through the charge dependent 
terms explicitly included in the two-body potential. The calculations have been performed using 
two different methods: the solution of the Faddeev equations in momentum space and the expansion 
on the correlated hyperspherical harmonic basis. The results are in agreement within 0.1% and can 
be used as benchmark tests. Results for the CD-Bonn interaction are also presented. It is shown 
that the and ^He binding energy difference can be predicted model independently. 


In the last years great efforts have been made to im¬ 
prove the description of the nucleon-nucleon (NN) inter¬ 
action. A new generation of potentials including explic¬ 
itly charge independence and charge symmetry break¬ 
ing (CIBjCSB) terms appeared. These interactions de¬ 
scribe the NN scattering data below Tiab = 300 MeV 
with a nearly perfect y^/datumRi 1. The CD-Bonn 
and Argonne ?;i8 (AVIS) [H interactions also provide a 
neutron-neutron (nn) force, which has been adjusted to 
the experimental nn scattering length, whereas the Ni¬ 
jmegen interactions are fitted only to proton-proton 
and proton-neutron data. Recently, the CD-Bonn poten¬ 
tial has been updated to CD-Bonn 2000 0 . In this paper 
we only present results for the AVIS and CD-Bonn 2000 
interactions. Both are quite different from each other in 
their functional form, but their description of the NN 
data is almost equally accurate. Therefore a comparison 
of the results will give insights into the model dependence 
of our understanding of the three-nucleon (3V) bound 
states. 

Following for example the notation of Ref. , all these 
NN potentials can be put in the general form 

viNN) = v^^{NN) + v^iNN) + v^iNN) . (1) 

The short range part v^{NN) of all of these interac¬ 
tions includes a certain number of parameters (around 
40), which are determined by a fitting procedure to the 
NN scattering data and the deuteron binding energy 
(BE), whereas the long range part is represented by the 
one-pion-exchange potential (NN) and an electromag¬ 
netic part {NN). 

For AVIS, v^^{pp) consists of the one- and two- 
photon Coulomb terms plus the Darwin-Foldy term, vac¬ 
uum polarization and magnetic moment interactions. 
The v^^{np) interaction includes a Coulomb term due 
to the neutron charge distribution in addition to the mag¬ 
netic moment interaction. Finally, v^^{nn) is given by 
the magnetic moment interaction only. All these terms 
take into account the finite size of the nucleon charge dis¬ 
tributions. The v^^{NN) for CD-Bonn is much simpler: 


v^^{pp) is given by the Coulomb force of point protons, 
whereas {np) = {nn) = 0. 

As it is well known, when these interactions are used to 
describe the ZN bound state, an underbinding of about 
0.5 MeV to 0.9 MeV depending on the model is ob¬ 
tained (see for example Ref.|||]). The local potentials 
lead to less binding than the non-local ones, a character¬ 
istic related to the bigger D-state probability predicted 
for the deuteron. Hence, it seems to be not possible 
to describe the A > 2 systems without the inclusion 
of three-nucleon interaction (TNI) terms in the nuclear 
Hamiltonian. Several TNI models have been studied in 
the literature mostly based on the exchange of two pions 
with an intermediate A excitation (for a recent review see 
Ref. [^). These interactions include a certain number of 
parameters not completely determined by theory, there¬ 
fore some of them can be used to reproduce, for example, 
the triton BE. 

In the following we show BE results for ^H and ^He. 
The AVIS interaction in conjunction with the Urbana 
IX (UIX) TNI [0 has recently been discussed by the 
Argonne group in a study of several bound states with 
3 < A < 8 [^. Here we perform a detailed calculation of 
the A = 3 system including total isospin states T = 1/2 
and 3/2. In addition, particular attention will be given 
to the BE difference D = i?(^H) — i?(^He) as a test of the 
CSB terms present in the interaction. The experimental 
value of this quantity is 764 keV, from which only 85% 
correspond to the standard Coulomb potential [^, |^, 0 - 
The remaining 15% should come from other CSB terms. 
A previous analysis of the contributions to D has been 
performed a decade ago j^, before the construction 
of the new interactions, which include for the first time 
the CSB terms in the fit to the NN data. Therefore a 
reanalysis is in order and might remove uncertainties due 
to an inaccurate description of the NN data. 

Because the Coulomb energy scales with the BE of ^H 
1^, we need ZN Hamilonians, which predict this observ¬ 
able accurately. This can be achieved with properly ad¬ 
justed TNI’s. Then the calculation of D requires reliable 
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Hamiltonian 

|A| 

B 

T 

Ps' 

Pp 

Pd 

Pt^3/2 

AV18 (T = 1/2) 



— 

7.618 

46.714 

1.295 

0.066 

8.510 

— 

AV18 (r = 1/2, 

3/2) 


— 

7.624 

46.727 

1.293 

0.066 

8.510 

0.0025 

AV18-bUIX(r = 

1/2) 


— 

8.474 

51.262 

1.055 

0.135 

9.301 

— 

AV18-bUIX(T = 

1/2, 

3/2) 

— 

8.479 

51.275 

1.054 

0.135 

9.301 

0.0025 

AV18 (r = 1/2) 



7.622 

7.616 

46.73 

1.290 

0.066 

8.510 

— 

AV18 (T = 1/2, 

3/2) 


7.621 

7.621 

46.73 

1.291 

0.066 

8.510 

0.0025 

AV18-bUIX(T = 

1/2) 


8.477 

8.470 

51.28 

1.051 

0.135 

9.302 

— 

AV18-bUIX(T = 

1/2, 

3/2) 

8.476 

8.476 

51.28 

1.052 

0.135 

9.302 

0.0025 


TABLE I: BE B, mean value of the kinetic energy T, S'-, P- and D-probabilities and the probability of the T = 3/2 state. 

The Pisa results are displayed in the first four rows. The last four rows show the Bochum results, in this case the modulus of 
the eigenvalue E is also given. All energies are given in MeV. The probabilities are given in %. 


Hamiltonian 

\E\ 

B 

T 

Ps' 

Pp 

Pd 

II 

w 

to 

AV18 (r = 1/2) 



— 

6.917 

45.669 

1.531 

0.064 

8.468 

— 

AV18 (r = 1/2, 

3/2) 


— 

6.925 

45.685 

1.530 

0.065 

8.467 

0.0080 

AV18-bUIX(r = 

1/2) 


— 

7.742 

50.194 

1.242 

0.131 

9.249 

— 

AV18-bUIX(r = 

1/2, 

3/2) 

— 

7.750 

50.211 

1.242 

0.132 

9.248 

0.0075 

AV18 (T = 1/2) 



6.936 

6.915 

45.70 

1.515 

0.065 

8.465 

— 

AV18 (T = 1/2, 

3/2) 


6.923 

6.923 

45.68 

1.524 

0.065 

8.466 

0.0081 

AV18-bUIX(T = 

1/2) 


7.759 

7.738 

50.23 

1.229 

0.132 

9.248 

— 

AV18-bUIX(T = 

1/2, 

3/2) 

7.746 

7.746 

50.21 

1.235 

0.132 

9.248 

0.0075 


TABLE II: Same as Table ^ for ®He. 


Hamiltonian 


®H 

3 

He 

AV18 

6 

keV 

-6 

keV 

AV18-bUIX 

7 

keV 

-7 

keV 

AV18 

6 

keV 

-6 

keV 

AV18-bUIX 

7 

keV 

-7 

keV 


TABLE III: Contribution of the proton and neutron mass dif¬ 
ference to the and ®He BE. The Pisa results are displayed 
in the Erst two rows. The last two rows show the Bochum 
results. 


solutions of the 3fV Schrodinger equation including these 
TNFs. For this reason we performed, as a by-product 
of our analysis, a benchmark calculation for the bound 
states of and ^He using the AV18-I-UIX potential 
model. The Bochum group solves the Faddeev equation 
in momentum space [|], whereas the Pisa group uses a 
decomposition of the wave function in pair-correlated hy- 
perspherical basis functions ©HI- Both methods were 
used to find BE’s to an accuracy of 2 keV, which means 
an accuracy better than 0.1%. Such a level of accuracy is 
nowadays routinely achieved for the 37V system W sev¬ 
eral methods using only NN interactions |||, |lj, |l^ 
Here we show that the same level of accuracy is obtained, 
when TNI terms are taken into account. 

We start considering ^H. The calculations have been 
done for three equal fermions using the isospin formal¬ 
ism. We used an averaged nucleon mass M with the 
value Ti^ /M = 41.471 MeV fm^ (the contribution of the 


Interaction term 

D 

nuclear CSB 

65 keV 

point Coulomb 

677 keV 

full Coulomb 

648 keV 

magnetic moment 

17 keV 

n-p mass difference 

14 keV 

total (theory) 

744 keV 

experiment 

764 keV 


TABLE IV: Contributions of the various terms of the inter¬ 
action to the ®H-®He mass difference D. The AV18-I-UIX po¬ 
tential has been used 



®H 

®He 



\E\ T 

\E\ T 

D 

CD-Bonn 2000 

8.005 37.64 

7.274 36.81 

0.731 

CD-Bonn 2000-fTM 

8.482 39.39 

7.732 38.54 

0.750 

Exp. 

8.482 — 

7.718 — 

0.764 


TABLE V: ZN BE’s \E\ for CD-Bonn 2000 with and without 
TM-TNI compared to the experimental values. Results are 
shown for ®H,®He and their BE difference D. Additionally, 
we show the kinetic energies T. All results are given in MeV 


n-p mass difference will be given separately). The AV18 
and AV18-I-UIX have been used to calculate the BE B, 
the mean value of the kinetic energy T as well as the 
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^|-E| D ann 

0.9990 8.474 742 -18.75 
0.9995 8.478 746 -18.86 
1.0000 8.482 750 -18.97 
1.0005 8.486 754 -19.08 
1.0010 8.491 759 -19.19 
1.0020 8.499 767 -19.42 

TABLE VI: Strength factor A for the nn force, resulting 
BE \E\ in MeV, the BE difference of ®He and D in 
keV and nn scattering length ann in fm. The calculations 
are based on the CD-Bonn 2000 potential modified by the 
strength factor in the ^So partial wave and the TM-TNI. 



a„„ [fm] 


FIG. 1: Difference of the ®He and BE’s D dependent on 
the nn scattering lengths a„ri. The crosses are based on the 
calculations shown in Table and the solid line is a linear 
fit to the crosses. The dashed, the dotted and the dashed 
dotted lines mark pairs of D and ann, which belong to the 
experimental D {ann is an estimation based on the linear fit 
in this case), the predictions of CD-Bonn 2000 and AV18, 
respectively. 


S', P- and £)-state probabilities. The results are given 
in Table || corresponding to two different calculations: 
i) for total isospin limited to T = 1/2 and ii) includ¬ 
ing also T = 3/2. The occupation probability Pt=zI 2 
of this state is given in the last column of Table The 
first four rows of the table show the Pisa group results, 
whereas the last four show the Bochum group results. In 
the later case, the BE \E\ is determined from the eigen¬ 
value spectrum of the Faddeev equations. Additionally, 
we present the absolute value of the expectation value 
of the Hamiltonian B. One sees that both values agree 
within the numerical accuracy. In Table the same set 
of results are given for ^He. For the T = 3/2 calcula¬ 
tions we find good agreement for the BE results and the 
wave function properties for both nuclei. The BE’s are 
in agreement within 2 keV or 0.05 %. The deviations 
for the wave function properties, especially for Ps', are 
slightly bigger, but remain below 0.4 %. This is below 


our numerical error bounds and confirms the reliability 
of both methods, even in presence of a TNI. The tables 
also reveal a small, but appreciable, contribution of the 
T = 3/2 state to the BE. Its inclusion produces 5 keV (8 
keV) more binding in (^He). It should be noted that 
the T = 1/2 results depend on the numerical method. 
The truncation of the Hilbert space to T = 1/2 leads to 
average pp {nn) and np matrix elements in the isospin 
t = 1 NN channels. This averaging is performed for 
the potential matrix elements in case of the Pisa calcula¬ 
tions, but for the t-matrix in case of the Bochum scheme. 
The small, but visible differences show that benchmarks 
to this accuracy require the comparison of fully charge- 
dependent calculations. This inconsistency also shows up 
in a visible deviation of \E\ and B for the T =ll2 Fad¬ 
deev calculations, because B is based on matrix elements 
of the potential wheareas |E| is based on the t-matrix. 

The contribution of the n-p mass difference is visible, 
but sufficiently small to be treated perturbatively. There¬ 
fore we show only perturbative estimates in Table HI. 
The positive sign in the tritium case indicates a slightly 
more bound system, conversely the ^He results slightly 
less bound. Again we find an encouraging agreement be¬ 
tween the Pisa and Bochum results. 


Taking into account the contribution of the n-p mass 
difference and averaging the Pisa and Bochum results, 
the final values of the BE’s for the AV18-I-UIX are: 
H(3 h) = 8.485(2) MeV and ^(^He) = 7.741(2) MeV. 
This is to be compared to the experimental values: 
Hexp(^H) = 8.482 MeV and Be,,p{^Be) = 7.718 MeV. 
Therefore, the AV18-I-UIX potential overbinds the tri¬ 
tium only by 3 keV, whereas the ^He is overbound by 
23 keV. This can be better analyzed looking at the pre¬ 
dicted BE difference D = 744 keV, which is 20 keV 
smaller than the experimental value. 

The calculations have been performed using the static 
Coulomb potential, i.e. the momentum dependence 
through a' has been disregarded ||^, so that our results 
correspond to o' = a. It has been estimated that the 
momentum dependence of a might contribute about half 
of the 20 keV shortage found in the value of D 
Therefore our results show that there is room from other 
contributions as for example CSB of TNI terms and a 
refinement of relativistic effects on D. 

The contributions to D of different parts of the in¬ 
teraction have been studied calculating the ^H and ^He 
BE’s omitting these parts and comparing to the full cal¬ 
culations. Note that this is not perturbative. The re¬ 
sults for the AV18-I-UIX potential including isospin states 
T = 1/2 and 3/2 states are collected in Table 0- We dis¬ 
tinguish i) the nuclear CSB terms, ii) the point Coulomb 
interaction. Hi) the complete pp and np Coulomb interac¬ 
tion, which includes the finite size charge distributions, 
the one- and two-photon terms and the Darwin-Foldy 
and vacuum polarization interactions, iv) the magnetic 
moment interaction, v) the n-p mass difference. 

One sees that the electromagnetic interaction on top 
of point Coulomb visibly contributes to D and cannot be 
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neglected. This raises the interesting question, whether 
CD-Bonn 2000, coming without an elaborate electromag¬ 
netic force, can also describe D. 

To this aim we performed 3N BE calculations using 
the CD-Bonn 2000 interaction. The results are given in 
Table Again, the NN interaction underbinds the 3N 
nuclei. Therefore we augmented the Hamiltonian by the 
Tucson-Melbourne (TM) TNI j^, |2^, |2^ . The strength 
of the original model has been adjusted to reproduce the 
experimental BE as described in |Q. It results the 
ttNN cut-off value A = 4.795 m,r- Again ^He is over¬ 
bound. The mass difference D, of 750 keV is slightly 
improved compared to our result for AV18 and UIX. 

In Ref. |2^ it has been observed that D is only sensi¬ 
tive to CSB in the S-wave. Unfortunately, nn scattering 
is only poorly known experimentally. There is only one 
datum for the scattering length, which is still controver- 
sal 1^, ^ 1^, ^ AV18 is adjusted to 

Arm = —18.82 fm, whereas for CD-Bonn 2000 a„„ results 
in —18.97 fm. Its charge-dependence is based on theoret¬ 
ical predictions of the full Bonn model Q. To pin down 
the origin of the difference of the predictions of both mod¬ 
els, we modified the ^Sq nn interaction of CD-Bonn 2000 
by a factor A and calculated the resulting nn scattering 
length a„ra, the BE and D. We found a strong linear 
correlation of and D shown in Table 0 and Fig. |. 

Moreover, the prediction of AV18-I-UIX perfectly fit into 
the results obtained from CD-Bonn 2000 and TM (see the 
dashed-dotted marks in the hgure). This shows that the 
dependence of D on the interaction can be traced back 


to different predictions for the nn scattering length. The 
very different treatment of electromagnetic interactions 
and the differences of the CSB in higher partial waves do 
not affect appreciable D. The deviation of the A^7V-|-TNI 
force prediction of D to the experimental value might be 
caused, as stated before, by CSB TNI terms not consid¬ 
ered in the present description or by relativistic effects. 
As a consequence we cannot estimate ann from our cal¬ 
culations of D. However, we would like to note that the 
a„„ —16.3 fm found in Refs. ^ would worsen our 

description of the 3N BE difference signihcantly. 

In summary we have calculated the ^H and ^He BE’s 
based on modern NN interaction models including TNI 
terms using two different numerical methods. Our re¬ 
sults showed the stability and reliablity of both schemes. 
Using only NN forces the BE’s are too small, calling for 
TNI terms. These lead by construction to the experi¬ 
mental ^H BE. We found that the BE difference of ^H 
and ^He is predicted nearly model independently. We 
could trace back the remaining model sensitivity to the 
differences in the predictions for the nn scattering length. 
However, uncertainties arising from CSB TNI terms and 
relativity do not allow us to extract the scattering length 
from the ^H and ^He BE difference. The model depen¬ 
dence arising from the different electromagnetic parts of 
the interactions is extremely small. 

A.N. acknowledge partial support from NSF grant# 
PHY0070858. Parts of the numerical calculations have 
been performed on the Cray T3E of the NIC in Jiilich, 
Germany. 
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